Ongoing

TO-DO

• Add equivalent of Raster.from_array() to RasterAccessor class (or RasterBase class?) and individual setting operations (for transform, crs, nodata, and area_or_point) to make all methods (reproject, etc) naturally work on both Raster and ds.rst,
• Ensure dual support of masked-array & NaN arrays for all methods (_reproject,, etc),
• Add delayed function (in _reproject, _interp_points and _subsample) by detecting automatically if input array is delayed.
• Add tests specific to RasterBase,
• Add tests specific to RasterAccessor functionalities (comparing to ds.rio).

Other Dask support to add (will be moved as issues for later PRs)

• The reduce_points function can copy the same logic as interp_points,
• The crop function using isel of Rioxarray,
• Potentially look at geocube for rasterize/polygonize support,
• But the proximity function would be a bit of work...

Resolves #383
Resolves #567